Matrix-element corrections to gg/qq 



Higgs in HERWIG 



G. Corcella^'"^ and S. Morettfi 

^CERN-TH, Switzerland; ^MPI fur Physik, Munchen, Germany; ^Southampton University, UK 
Abstract 

We describe the HERWIG implementation of real matrix-element corrections 
to direct Higgs hadroproduction at Tevatron and Large Hadron Collider (LHC) 

. . . and compare it to other approaches existing in literature and describing the 

^ . transverse momentum distribution of the Higgs boson. 

o 

<N ; 1. THE HIGGS TRANSVERSE MOMENTUM 

■ In order to investigate Higgs boson production via gg Higgs (see Ref. d), one needs to account 

pL| \ for multi-parton radiation for the sake of performing trustworthy phenomenological analyses |2l|3l|4l. 

ff^ \ Standard Monte Carlo (MC) algorithms EJ-Q describe parton radiation in the soft and/or collinear 

approximation of the parton shower (PS), but can have regions of phase space, so-called 'dead zones', 
where no radiation is allowed. Here, one can however rely on higher-order tree-level results, as in 
^ \ this region the radiation is neither softly nor coUinearly enhanced. Several methods have been recently 

^ ' suggested in order to match PS and fixed-order matrix elements (MEs) ||8||9l, also including the virtual 

one-loop terms fT0l - fT2l . 



2. THE HERWIG IMPLEMENTATION 



(N 
O 

^ . In this note, we briefly mention that the same strategy which has already been used to implement real 

\ ME corrections to e+e~ annihilation into quark pairs fljll . Deep Inelastic Scattering (DIS) lll4l . top 

I ■ quark decay lITSl and vector boson hadroproduction (T^ has now also been adopted for the case of 

Higgs hadroproduction via gluon-gluon fusion and quark-antiquai^k annihilation ^V.], in the context of 
(— I \ the HERWIG event generator |5l|6l. That is, the dead zone is here populated by using the exact next-to- 

leading order (NLO) tree-level ME result and the PS in the already-populated region is corrected using 
, , the exact amplitude any time an emission is capable of being the hardest so far. 

X 

. p. , 3. NUMERICAL RESULTS AND COMPARISONS 

The MEs squared for the real corrections to gg ^ H that we have used can be found in (TE\ . where top 
mass effects are fully included. The real NLO corrections to — > are instead rather straightforward: 
the formulae we used can be read from Eq. (3.62) of fTOll with appropriate Yukawa couplings and cross- 
ing. In the new HERWIG default version, in line with |[T6l . ME corrections use the Higgs transverse 
mass = gf^ + mj^ as the scale for as and for the Parton Distribution Functions (PDFs) while the 
99i QQ ^ H contributions use rnjj. We shall also assume that the intrinsic transverse momentum of the 
initial-state partons is equal to qT;mt = 0, the HERWIG default value. 

By adopting the HERWIG defaults, we first consider Higgs production at the Tevatron and the 
LHC within the MC itself, by plotting the qx distribution with (solid histogram) and without (dotted) 
ME corrections: see Fig. ^ Beyond qt ^ rriu/'i the ME-corrected version allows for many more 
events. In fact, one can prove that, within the standai^d algorithm, qt is constrained to be qT < rriH- 
At small qx the prediction which includes ME corrections displays a suppression. By default, after the 
latter are put in place, the total normalization still equals the LO rates. Hence, it is obvious that the 
enhancement at large qr impUes a reduction of the number of events which are generated at small qx 
values. 

In Fig. 121 (left plot) we present the improved HERWIG spectrum (solid) for the LHC, along with 
the result obtained running the so-called 'i7 -i- jets' process (dotted), where the hard process is always one 
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Fig. 1: Higgs transverse momentum distribution according to HERWIG with (solid) and without (dotted) ME 
corrections, at Tevatron (left, ^/s^^ ~ 2 TeV) and LHC (right, VSpp = 14 TeV). We have set the Higgs mass to 
■niH = 115 GeV. 
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Fig. 2: Left: comparison of ME-corrected HERWIG predictions (solid) to the 'iJ + jets' result from 1181 (dotted). 
Centre: comparison of ME-corrected HERWIG predictions (solid) to the NLO and resummed calculation of 1201 
(dotted). Right: comparison of ME-corrected HERWIG predictions (solid) to the MC@NLO results from the code 
described in Ref. II21I (dotted). Here, qq ^ H processes have been turned off. 



of the corrections to gg — > H. In order to perform such a comparison, we have turned the qq ^ H hard 
process off, as 'H + jets' in HERWIG does not cuiTcntly implement the corrections to quark-antiquark 
annihilation. Furthermore, we have chosen grmin = 30 GeV for the 'H -i- jets' generation. As expected, 
at small qr the two predictions are fairly different but at large transverse momentum they agree well. 

In Fig. |2l (centre plot) we compare the new HERWIG version with the resummed calculation 
of Ref. |20|. For the sake of comparison with HERWIG, which includes leading logarithms and only 
some subleading terms, we use the results of ||20| in the NLL approximation (rather than the default 
NNLL one), matched to the NLO prediction. In order for such a comparison to be trustworthy, we 
have to make parameter choices similar to lEOl : namely, we adopt a top quark with infinite mass in 
the loop and inn = 125 GeV, with as and PDFs (both from HERWIG defaults) evaluated at m^. 
While the normalization (LO in HERWIG, NLO in Ref. iGOl) and the small-g^ behaviour of the two 
curves are clearly different, the large transverse momentum predictions are in good agreement, as in both 
approaches it is the real NLO ME that dominates the event generation at large qx- 

Finally, in Fig. |2l (right plot), we compai^e the results of standard HERWIG after ME corrections 
with the so-called 'MC@NLO' event generator (version 2.2) of Ref. ll2Tl . the latter implementing both 
real and virtual corrections to the hard-scattering process, in such a way that predicted observables (in- 
cluding normalization) are correct to NLO accuracy. As version 2.2 of the MC@NLO includes only the 
corrections to Higgs production in the gluon-fusion channel, we again have turned the quark-annihilation 
process off in our routines. As observed in the comparison with the resummed calculation, the two spec- 
tra differ in normalization and at small q^, but agree in the large-transverse-momentum region. 



4. CONCLUSIONS 

Between the described implementation and the one available within the MC@NLO option, we believe 
that HERWIG is presently a reliable event generator for (direct) Higgs production from parton fusion 
at hadron colliders both at small and large transverse momentum. In fact, all currently available ME 
corrections will play an important role to perform any analysis on Higgs searches at present and future 
colliders. In particular, the option described here may be the most convenient choice for when the phase 
space is Umited to transverse momentum values such that qt ^ mn- 
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